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ABSTRACT 


This work finds out the stresses developed in the critical sections of the two models of the three 
wheelers Vikram chassis frame at selected cross sections. Computer modeling of the two chassis 
has been done. Using this computer model finite element analysis of the two chassis has been 
undertaken to determine and compare the stresses developed in them at critical sections. Analysis 
has been done for three different dynamic loading conditions. For one loading condition 
experimental stress analysis at some selected points has been performed to verify the FEA model 
and results for both the models. 
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CHAPTER - 1 


INTRODUCTION 


1.1 PROBLEM DEFINITION 

Three wheelers are popular city transport vehicles specially for short hauls of passengers and 
loads. One such three wheeler “VTKRAM”' is manufactured by Scooters India Limited, 

Lucknow. This three wheeler has a major share of market in U.P.. 

Vikram was initially imported from an Italian company - Innocenti. This vehicle was designed 
with a two-stroke petrol engine. To make it suitable to the Indian requirements, the manufacturer 
m India has modiRed it by (a) replacing petrol engine with a diesel engine and (b) by increasing 
the length of the chassis. 

The initial two stroke petrol engine of the three wheeler Vikram was very small in size, half the 
weight of the diesel engine now being used and it ran at much higher engine speeds as compared 
to the present one. Being a two stroke engine it had lesser unbalanced forces. The present engine 
is single cylinder four-stroke diesel engine primarily designed for stationary agricultural needs 
such as pumping water, running the thrashers etc. Because the engine has been designed to be 
anchored at one place and not for a vehicle, the unbalance forces are considerable along with the 
high noise level. Added to it are the bad road conditions, rash dnving, and overloading due to 
which large forces are transferred through the road to the chassis frame. 
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Fig. 1.1 Top view of three wheeler Vikram (Side mounting) 


Fig. 1.2 Front view of three wheeler Vilcram (Side mounting) 
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Fig. 1.3 Top view of three wheeler VDcrarn (Central mounting) 


Fig. 1.4 Front view of three wheeler Vikram (Central mounting) 
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Fig. 1.5 Three wheeler Vikram diesel engine 
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Due to dynamic forces and static loads, failure has been reported at three different sections of the 
chassis during the guaranty period of the vehicle. To over come this problem, the manufacturer 
has redesigned the old chassis. In the new design, the engine has been located in the middle of the 
chassis instead of the side. With the changed engine location, they have also changed the load 
carrying members. Since these changes have been incorporated very recently the manufacturer 
has no feed back concerning the new design. 

Thus, the problem concerns the following: 

(a) Computer modeling of the two chassis (with the old and new designs) acted upon by dynamic 
and static loads. 

(b) Finite Element Analysis (FEA) of the two chassis models to determine and compare stresses at 
critical sections under the loading conditions. 

(c) Experimental stress analysis at some critical sections verifying the FEA model and results. 

1.2 LITERATURE SURVEY 


In the moving automobile, different loads act on its parts and units. The magnitude and nature of 
these loads depends on wheel-road interaction, driver’s action on controls, engine operating 
conditions and other factors. The reliability of an automobile depends on actual loads that act on 
it when in use. The automobile is the most complex object for determining load conditions, as its 
service conditions are very diverse and vary with time. Variable factors combine differently; the 
magnitude of stresses and the nature of their distribution in different parts vary within wide limits 
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The forces causing these stresses may be constant (such as weight, forces which appear when the 
parts are tightened at the time of assembly, and the like) or variable (starting and braking forces, 
forces caused by manufacturing inaccuracies, road and operating conditions, resistance forces, 
dynamic forces, and others). Two types of strengths are recognized; static and fatigue strength. 
Static strength is the ability of a part to resist failure under the action of short time maximum 
loads. Fatigue strength is the ability of a part to resist failure under the action of repeatedly 
occurring loads, continuously varying or cumulative. 

Since long, designers have been trying to design the automobile to make it for smoother ride, 
more comfortable, safer and stronger. Costin and Phipps [ 2 ] have tried to do chassis frame’s 
stress calculation using force equilibrium method. They have found out the stresses developed in 
the different loads carrying members of their racing car for the static case. 

Ali and Hedge [ 1 ] have done a FEM idealization of truck chassis frame and used it to predict 
deflections under various types of loads. They have also built a test rig to support the chassis 
frame to find out the deflections of the frame relative to a specially supported reference frame for 
out of plane bending, torsion and lateral loadings. 

Sharman [ 8 ] has performed an experiment in which he recorded the strains at critical points on 
the chassis of a 16 ton GVW semi-trailer, operating on various road surfaces. He has used the 
statistical analysis of the records in terms of design factors and load spectra to determine the 
maximum static loads and have estimated the fatigue life. 
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Generally, the automobile industry have made various efforts for the better design of a chassis 
frame, no such data is available with the manufacturer of the particular three wheeler Vikram. 

s. 

The only work available on the three wheeler Vikram is done by Joshi [ 3 ]. For three wheeler 
Vikram’s Side mounting frame, he has tried to bring out an improved design of the engine- 
mounting system which encounters least number of resonating phases during operation and thus 
looks for a design that results in less breakage of the exhaust pipe. 


1.3 THE PRESENT WORK 


The present work aims at predicting the behavior of the strength of the critical sections at the 
selected cross sections (node points) for both the models by doing finite element analysis for the 
following three different cases: 

CASE I: When the engine is idling and vehicle is stationary. 

CASE II: When the vehicle is running at 40 km ./hr. in the fourth gear and its front wheel goes 
into a ditch of shape of half sine wave of 2 m wave length and 150 mm amplitude. 

CASE HI: When the vehicle is running at 40 km./hr. in the fourth gear and its both the rear 
wheels go into a ditch of shape of half sine wave of 2ra wave length and 150 mm amplitude. 

This work gives time Vs factor of safety plots for these selected points using maximum shearing 
stress theory. It compares two models at the critical sections for their strengths for all the three 



cases. It also verifies the finite element model and results by experimental analysis on the selected 
points, for case I, for both the models. 

Chapter 2 deals with the Finite element method formulation and those features of NIS A related 
to the present work. The geometric modeling of the two types of chassis frames has been done in 
chapter 3 for pre-processing with NIS A. The analysis results from the modeling are stated and 
discussed in Chapter 4. Experimental details and observations for the chassis frames are given in 
chapter 5. Finally, chapter 6 recommends the scope for further work that can be undertaken on 
the vehicle. 
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CHAPTER - 2 


NISA - AN OVERVIEW 


2.1 FINITE ELEMENT METHOD - A DEFINITION 

Finite element method [ 10 ]is a powerful numerical technique for approximate solution of 
continuum mechanics and can be regarded as an extension of the Rayleigh-Ritz Procedure where 
piece-wise continuous trial functions are used. The method is approximate, since a continuum 
with infinite numbers of degrees of freedom is replaced with a discrete system with finite number 
of degrees of freedom. The method involves subdividing a continuum into a finite number of 
regions called elements connected at finite numbers of points called nodes to which loads are 
applied. An approximate admissible solution is constructed over the assemblage of elements, and 
the solution continuity is maintained at the inter-element boundaries. 


2.2 NISA 


NlSA(Numerically Integrated elements for System Analysis) is a general purpose finite element 
program developed and maintained by EMRC to analyze a wide spectrum of problems 
encountered in engineering mechanics. 
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2.2.1 GOVERNING EQUATION OF MOTION 


The equations of motion describing the dynamic behavior of a structural system may be obtained 
from the extended Hamilton’s principle for elastodynamics as; 


5j(t/-T)dr-j5Wt = 0 

fo ^0 


( 2 . 1 ) 


where 

5 = first variation 

to,tf= two arbitrary time points at which the first variation vanishes 
U= strain energy in terms of geometrically compatible displacements 
T= kinetic energy in terras of geometrically compatible displacements 
5 W= virtual work of the external forces acting on the system during virtual displacements 

The external forces may be conservative or nonconservative. A force is said to be conservative 
when its components are derivable from a potential function Nonconcervative forces, e.g. 
damping and follower forces, are deformation dependent, their magnitude and/or direction depend 
on the deformation. During the virtual displacements, the variations of the nonconservative forces 
are ignored. The virtual work 5 W, therefore ,can be written as: 
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5W = -§y + 5W„, (2.2) 

where V is the potential of the apphed conservative forces and 5W x is the virtual work of the 
nonconservative forces during virtual displacements. Substituting Eqn. (2.2) into Eqn.(2.1),one 
obtains 

5j(n-^^f-j5W„,Jr = 0 (2.3) 

^0 ^0 


where is the total potential energy of the system .Above equation reduces to the static case 
when the kinetic energy vanishes and the remaining terms are invariant with time. 

Here viscous damping forces will be the only nonconservative forces considered. Introducing the 
discretization into the kinetic energy and the virtual work of nonconservative forces: 


n=n(“)=in‘(““’) (2^) 

tf=l 


T = r(if) = X ) = ijpu''udv 




(2.5) 
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( 2 . 6 ) 


e=l 


where u and u are the nodal displacement and velocity vectors, respectively, and and 
are respectively the consistent mass and damping matrices for element e ,given by: 




(2.7) 


C(‘) ^ f 

Jv ^ 


( 2 . 8 ) 


where and c/'^ are the mass density and the viscous damping coefficient for element e, and N 
is the matrix of shape functions. The minus sign is used in Eqn (2.6) since the damping 
force(viscous damping which is proportional to the velocity is assumed) opposes the motion. 
Using Eqns.(2.4) through (2.6) in the variational principle (2.3),we obtain: 

J + ^]dt = 0 (2.9) 

■"o du du du 

where p is the vector of applied loads and fd is the damping force vector (assembly of fd^'^ of 
Eqn. 2.6) associated with the generalized coordinates u . Integrating the second term (involving 


13 



« ) of Eqn.(2.9) by parts and observing that 5« vanishes at t = to and t = tf and arbitrary in the 
open interval (to,tf), we obtain: 


d dT ^ dU 
dt du du 


dT 

du 


( 2 . 10 ) 


which is a statement of Lagrange’s equations. From this we can derive the governing differential 
equations of motion(the dynamic equilibrium equations) 


M 11 + Cu + K u = p(t) 


( 2 . 11 ) 


where M, C and K are respectively the global mass, damping and stiffness matrices(made up by a 
proper assembly of the element matrices), p(t) is the time-dependent apphed force vector, and u, 
u and u are the nodal acceleration, nodal velocity and nodal displacement vectors, 
respectively(the bar is dropped for convenience, i.e., u=m ). 


2.2.2 DIRECT TRANSIENT ANALYSIS 


As shown, the governing differential equation of motion can be written as: 


Mu + Cii + Ku = p(t) 


The corresponding initial conditions can be written as: 
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u(o) = do 


and 


u(o) = Vo 


where do and Vo are the given initial displacement and velocity vectors, respectively. Direct 
transient analysis is based on direct integration of the equation of motion, i.e. no prior 
transformation or eigenvalue analysis is required. 

In NIS A, Newmark method is used for direct transient analysis, which can be summarized by the 
following equations: 


Mu'n + CUn + KUn = Pa (2. 12) 

Un+I = Un + At Un + ( y - P)At^ u'a +■ ^At^ Un+1 (2. 13) 


Ua+l = Ua + (l-Y)At Un + YAtUn+i (2. 14) 


where Ua,Ua and Un are the known solutions for u, u and u at time step n, respectively. The 
parameters P and y determine the stability and accuracy characteristic of the algorithm. In NISA, 
the Newmark method is implemented with displacement form. The predictor are defined as: 
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• * ■f 

Un+1 = Un + (1 - Y)At Un 


Un+1 = Un + At Un + (3 - | 3 )At^ Un 


With this and Eqns. (2.13) and (2.14) the velocity and the acceleration of the current time step 
can be written as: 




(2.15) 


• ^ 


y o 


PAr 


u. 


y 


n+l 


PAr 


TI + l 


(2.16) 


Substituting these into Eqn. (2. 12), we have the following recursive equation to solve for 
displacement Un+i: 


(—^M + -^C + K)u^, = p_+(— ^M + -^ 
^pAr' pAr pAr" pAr 




n+l 


(2.17) 


By defining the effective stiffness matrix. 


a:* 


PAf" 


■M + 


-^C + K 
pAr 


and the effective load vector. 




pAr" 


pAr 


C)^ 


/Tv 


’/I + l 


Cu 


'/r+l 
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then, Eqn. (2.17) becomes 


K*Un+i=p* (2-18) 

which is in the fonn of the familiar matrix equation for static analysis that can be solved for the 
displacement vector Un+i. Eqns. (2.15) and (2.16) can then be used to solve for acceleration and 
velocity. 

2.2.3 BODY FORCE 


In NIS A gravity loads and inertia loads(due to angular velocity and Unear angular acceleration) 
may be appUed to any element which possess a non zero mass density. With reference to Fig. 2.2, 
the XYZ system is the rotating frame of reference which rotates with an angular velocity (O and 
angular acceleration a about the X'Y'Z' system. At the instant shown, the two systems are 
parallel to each other. A Unear acceleration field at is also applied to the two systems. 

The absolute acceleration of a generic point p can be written as: 

ap = at + u + AT + n£2r + 2nu 


where 


and 


r — tpo “ to'o 
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el: uniform linear 
tccelcrition field 



Fig. 2.2 Acceleration at a point due to linear and angular motion 
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at = the absolute acceleration of the moving frame of reference X'Y'Z' due to the linear 
acceleration field. 

u = the stretching-term acceleration (i.e., relative) as viewed by an observer attached to the 
X'Y'Z' system. 

Ar = the tangential acceleration due to the angular acceleration a. For steady motion about a fixed 
axis, this term goes to zero. 

nnr = the centripetal acceleration (due to the angular velocity of the moving axes.) 

2nu = the Corriolis acceleration(due to the interaction between the angular velocity and the 
stretching velocity). 

For linear static analysis, the terms involving u and u are not considered, and the generic point is 
considered on the undeformed configuration. 
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2.2.4 NISA ELEMENTS 


Following is the brief description of the NISA elements used for modeling: 

2.2.4. 1 3-D Beam Element (NKTP = 12) 

This element (Fig. 2.3) is a 2-node pnsmatic 3-D beam element. The formulation includes 
stretching, bending and torsion effects. The transverse shear deformation effect is included as an 
option. The beam vertices may be offset from the corresponding nodal points and the centroid 
may be offset from the shear center. The element has end release capability. The deformation is 
characterized by three translations (UX, UY, UZ) and three rotations (ROTX, ROTY, ROTZ). 
The local x-axis of the beam is along the centroidal axis. The local y and z axes are user defined 
and are not necessary principal axes of the cross section. 

2.2.4.2 3-D SPAR ELEMENT(NKTP = 14) 

This element (Fig. 2.4) is a 2-node uniaxial tension-compression 3-D line element and is used to 
model space trusses. The element may be oriented any where in space. The element has three 
translational degrees of freedom per node (UX, UY, UZ). 
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Fig. 2.3 Element configuration and orientation for 3D beam element, NKTP=12 
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NORDR = 1 



Fig. 2.6 Element configuration and face numbering convention for top and bottom srufaces of 

sheU element, NKTP=20, NORDR=l 


Y 


^ X 



Fig. 2.7 Element configuration for 3D concentrated mass element, NKIP=26 



2.2A.3 3-D TRANSLATIONAL SPRING ELEMENT(NKTP = 17) 


This element (Fig. 2.5) is a 2-node uniaxial tension-compression mass less spring in three 
dimension. The element may be oriented anywhere in space. The element has three translational 
degrees of freedom per node(UX, UY, UZ). The two nodes defining the spring may be 
coincident, in which case a direction vector is specified indicating the direction of the spring axis. 

2.2.4.4 3-D SHELL ELEMENT(NKTP = 20) 

The 3-D shell element (Fig. 2.6) includes membrane, bending and transverse shear deformation 
effects and is suited for modeling moderately thick to thin shell structures. The element has six 
degrees of freedom per node ( UX, UY, UZ, ROTX, ROTY, ROTZ), but it posses no rotational 
stiffness about the normal to the shell surface. 

The element can be shaped as a 4 to 12 node quadrilateral, or a 3 or 6 node triangle. 

2.2.4.5 3-D CONCENTRATED MASS ELEMENTS(NKTP = 26) 

The 3-D concentrated mass element (Fig. 2.7) provides a mean to represent the inertia properties 
of the model and to apply external loads due to body forces(e.g. gravity). The inertia properties 
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may also be represented by the mass density of the structural elements. This element has three 
degrees of freedom(UX, UY, UZ). 

2.2.5 INPUT DATA SETUP FOR LINEAR DIRECT TRANSIENT DYNAMIC 
ANALYSIS 


Fig. 2.8 shows a global description for the input data setup for a typical linear direct transient 
dynamic analysis run. The first data block (executive commands) must flag the analysis type as 
linear direct transient and should include appropriate or applicable commands for file saving, 
element resequencing etc.. 

The second data block (the model data) should start with the *TTTLE data group in order to 
allow printing of the problem as a heading on all pages of the output file. The model data then 
describes model characteristics, as for example: element definition (*ELEMENT), nodal 
coordinates (*NODES), material properties(*MATERIAL),etc.. Some kinematic constraints are 
considered to be model characteristics and properties in this block, namely: the rigid element 
data(*RIGLINK), the multi-point constraint equations(*MPCEQN) and the couple displacement 
data (*CPDISP). Other general groups such as definition of vectors and sets (* VECTORS and 
*SETS, respectively) are also defined in this block. 
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The third data block (the analysis data) should start with *EVENT data group. An event tide is 
available on the *LCTi'iLE data group. The block then describes various loading, kinematic 
boundary conditions and initial conditions of the structure. The initial condition data (*INTTAL) 
can appear before the *EVENT data group or in the first event. Loading include, for example, 
concentrated force (*CFORCE) and body force(*TBODYFORCE). Kinematic boundary 
conditions only contain the specified displacement data (*SPDISP). The block defines various 
output and print control options. It should be noted that the specified displacement data, 
*SPDISP, cannot be changed from one event to another. 

The last card of the input deck must always be the input data terminator, the *ENDDATA data 
group. 
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Fig. 2.8 Input data setup for linear transient dynamic analysis 
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CHAPTER - 3 


MODELING OF CHASSIS FRAMES FORNISA 


3.1 INTRODUCTION TO THE CHASSIS FRAMES 


Both the present chassis structures can be looked upon as a space frame consisting of 3-D beam 
elements, 3-D spar element, 3-D translational spring element, 3-D general shell element and 3-D 
concentrated mass element Location of joints and directions of members are completely arbitrary, 
and the structure is designed to carry forces that are general in space. 


3.1.1 SIDE MOUNTING MODEL 


Fig. 3.1 shows the chassis of the three wheeler ‘VUCRAM’ (side mounting). It is made of a 
central tube through out the length of the vehicle. This tube can be called the central main element 
of the entire vehicle. Two channels running parallel to the tube are mounted to the tube through 
cross member, which is perpendicular to the tube attached to the tube. These cross members are 
also channel section having the same cross sectional area as that of side channel. A platform made 
of same cross section channel is supported by these side channels as well as two leaf springs in the 
rear part of the chassis(one on the each side). Other end of the leaf spring is attached to the axle 
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joining the two reax wheels. Engine and gear box is mounted between central tube and one of the 
side channels, on plate at six different supporting points. Front tire of the vehicle is attached to 
spring which is connected to a semicircular fork through spar elements. The middle point of this 
fork is attached to the central tube through another tube. Cross sectional properties for this model 
are given in the tables 3.1 to 3.5 and in the Fig. 3.10. Section 1, section 2, section 3.1 and section 

3.2 (Fig. 3. 1) are critical sections for this model. 

3.1.2 CENTRAL MOUNTING MODEL 


Fig. 3.2 shows the another model of the chassis of the three wheeler ‘VIKRAM’( central 
mounting ). Front part of it is made of the central tube and is similar to the side mounting model. 
Central tube is cut when it becomes honzontal. It is joined to two tubes running parallel to each 
other in the horizontal plain, through cross members at two different points. These cross members 
are also tubes of same diameter. Rear platform is supported on parallel tubes as well as on two 
leaf springs (one on each side). Other end of the leaf spring is attached to the axle; joining the two 
rear wheels. Engine and gear box is mounted between two parallel tubes, on plates at six different 
supporting points. Cross sectional properties for this model are given m the tables 3.6 to 3.10 and 
in the Fig. 3.11. Section 1, section 2.1, section 2.2 and section 3.1 to section 3.4 (Fig. 3.2) are 
critical sections for this model. 
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3.2 ENGINE SPECIFICATIONS 


Following are some particulars of the VIKRAM engine 

Make: Greaves India Limited, Aurangabad 

Type: Single cylinder, 4-Stroke air cooled diesel engine 

Cylinder Bore: 85 mm 

Stroke Length: 90 mm 

Cubic Capacity: 510 cm^ 

Max. Output: 1 1 H.P. at 3000 rpm 
Max. Torque: 3 kg-m at 2000 rpm 
Engine weight with gear box assembaly: 102 Kg 
Speed range of the engine: 850-3200 rpm 
Engine Power Pack Dimensions: 

Engine height: 562 ram 

Engine length: 753.6 ram 

Engine width: 232 mm at mounting points 



3.3 OTHER MISCELLANEOUS SPECIFICATION 


Gear Ratios: 

Different gears ratios are given as bellow; 

T‘ gear = 1:5 
gear = 1: 2.93 
gear = 1: 1.84 
4* gear = 1: 1.29 
Reverse gear = 1; 5.34 
Secondary transmission (differential unit) 

Ratio = 1 : 4.1 

Wheel dia = 488 mm 

Material of the chassis structure: 

Entire chassis structure is made of oil quenched and tempered carbon steel with carbon contents 
up to 0.2%(max.). Various properties of the material are as follows: 

Tensile strength: 640 Mpa 
Yield stress: 440 Mpa 
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Specific Weight: 7.85 * 10"^ Kg/mm^ 
Poisson’s ratio: 0.292 
Modulus of rigidity: 0.79 * 10^ Kgf/cm^ 
Young’s modulus: 2.07 * 10® Kg/mm 

3.4 ENGINE FORCES AND MOMENTS 


At engine mount system following external forces are introduced due to the weight and moments 
associated with the engine. 

i. Weight of the engine: 

The engine mass including the mass of the gear box is 102 Kgs. As a result 17 Kgs load is 
introduced at the six engine mounting points. But the magnitude of this force is small as compared 
to the dynamic forces taken up by the chassis. 


ii. Reciprocating inertia forces, Free , in the engine: 


These forces are given [ 5 ] as. 
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Kec = V{cos(cor) + 


cos(2a)?) 

n 


} 


(3.1 


where, 

Mrec : reciprocating masses ( mass of piston = mass of gudgeon pin + 1/3 mass of con. rod) 
CO: angular velocity of the engine crank shaft 
r: radius of crank 

n: connecting rod length / crank radius 
A^: Area of piston 
Ip: connecting rod length 
r: crank radius 

n: Ip/r ( equal to 4 for diesel engine ) 

For Aluminum piston [ 5 ] 


0.43 ( 0.00415 D cylinder ) 
=1.095 Kg. 
h4coo>rod" 250 

= 1.42 Kg 

Mrec = 1.095 + (1.42/3) 

=1.568 Kg 
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This force in the present system acts vertically. In balanced engine, half of the magnitude of 
primary forces are balanced by placing the required counter weight on the crank shaft. But in 
doing this balancing, another force in horizontal direction is generated, given by, 

Fhori= (1/2) M^co^r sin(cDt) (3.2 

iii. Forces due to rotating masses, Frot j of the engine: 

It is assumed that the present engine is balanced to the extent that C.G. of the crank shaft and the 
crank lies along its axis. So now the rotating forces, will stem from the portion of the mass of the 
connecting rod, Mrot , which can be assumed placed at the big end, will be responsible for 
centrifugal forces. 


Mrot — (2/3)Mcon rod 
=0.94 Kg 


leading to. 


Frot = Mrot CD^ r sm(cot) 


(3.3) 




z 


Fig. 3.3 Schematic diagram of Engine - mount 
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Pislon-cylLnder arrangement 



[F l»n 9 



(b) 


(c) 


Forces coming on the engine Forces coming on the engine-mount 

Fig. 3.4 Forces resulting in the engine 
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Rc-:'king couple 



Fig. 3.5 Nature of rocking couple for 4-stroke engine at 1200 rpm 
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iv. Rocking Moment of the engine: 


In the engine, at the end of compression stroke, due to explosion inside the cylinder, a lateral 
force normal to the cylinder axis acts on the piston, besides the thrust in the connecting rod. As £ 
reaction to this lateral force , and opposing force acts on the cylinder walls. Figs. 3.4 b and 3.4 c 
shows the two forces F and F tan0 acting on the engine and the engine mounting supports 
respectively The force F acts upwardly on the cylinder head, is balanced by equal , down ward 
directed force on the crankshaft bearing. The summation of vertical forces for a engine structure 
is thus zero. The lateral force F tan0, spaced L distance apart, results in a moment about crank 
shaft axis. A.S a result , whole of the engine rocks. The nature of this moment is shown in 
Fig.3.5. This rocking moment is assumed to be composed of a constant component and a 
harmonic component. The engine being a 4 stroke engine, the frequency of harmonic component 
is half the engine speed, as this moment is dependent on the explosion frequency i.e. firing 
frequency. The rocking moment is assumed to be maximum at 20° ( assuming maximum pressure 
of about 60 bar inside the cylinder at around 20° [ 5 ]) crank angle as the explosion pressure is 
maximum around that crank angle. 

Maximum explosive force generated inside the engine, 

Fmax — Pmax A.p 

= ( 60 * 10^ ) 0.0056 

= 34344 N 
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Rocking moment. 


= (1/2) Icon tan(20°) [ 1 + sin( ox/2)] 

= (1/2) (34344) 0.18 tan(20°)[ 1 + sm(ox/2) ] 

= 1124.5 [ 1 + sin(cot/2)] (3.4 

So the net vertical force, 

Fy = (^M„, + M„,)o) Vcos((or) 

=1.724 co\cos(cot) (3.5; 

and the net horizontal force, 

K = iKot -|M„Jo)Vsin(cor) 

= 0.16 X sin(cDt) (3.6) 
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LliAF SPRING MODELING 


Two leaf springs support the rear platform of the chassis, one on each side. For the leaf springs 
the distance between two supporting points is 775 mm and the vertical distance of supporting 
pomt to the point which is attached to the axle is 270 mm. 

While raodehng for leaf spring it is assumed that both ends of the leaf spring are fixed to the 
platform. During modeling each leaf spring of stiffness K is replaced by two transitional springs o 
stiffness K each as shown in the Fig. 3.6 

From simple calculation we obtain the relation between K and K as: 

K’ = (K/2) sec^e 
if K =5.13* 10^ N/m 
then K = (5.13 * 10^2) (1+ 775/(2*270)) 

= 6.25 * 10^ N/m 
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3.6 BOUNDARY CONDITIONS 


3.6.1 DISPLACEMENT BOUNDARY CONDITIONS 

Since NISA does not provide the facility to model for tyre of an automobile the displacement 
boundary condition has been applied on the node at which tyres are mounted on the axle and the 
stiffness of the tube has not been considered. The following displacement boundary conditions 
have been applied for three different cases: 

Idling: 

CASE I 

In case I a stationary vehicle with idling engine is simulated. For this aU the six degrees of 
freedom of the rear tires’ mounting node are assigned zero. Since front tire is mounted on spar 
element its all three translational degrees of freedom at that node are also zero. 

Ditch: 

While simulating a ditch it has been assumed that when the vehicle goes into the ditch it is 
running at the linear speed of 40 km./hr. in the fourth gear. A ditch has been simulated with a 2 m 
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wavelength half sine wave of 15 cm amplitude. Wheel does not loose contact with the ditch 
surface while going through it 
Related calculations are: 

Speed of the vehicle = V = 40 Kms/hr = 11.11 m/s 
Time penod = 2/11.11= 0.18 sec 

Time taken by wheel while traveling through ditch = 0.09 sec 
Frequency = f = 1/T = 1/0.18 
.-. CO = 2 7C f = 34.9 rad/sec 
for motion through ditch 
z = 150sin(cot) 

CASE 11 

In case 2 a vehicle with front wheel passing over a ditch has been simulated. In this case front 
wheel’s mounting node has been given displacement of 
z = 150 sin(cot) 
where co = 34.9 rad./sec. 

in the z direction. Other two degrees of freedoms for it are zero. All the six degrees of freedom of 
rear wheel’s mounting node are also zero. 

CASEHI 

In this case rear wheels’ mounting nodes have been given displacement of 
z = 150sin(cot) 
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where co = 34.9 rad./sec. 


in the z direction. All the other degrees of freedom for it are zero. All the three degrees of 
freedom for front wheel’s mounting node are also zero. 

3.6.2 FORCE BOUNDARY CONDITIONS 

Static forces: 

3-D point mass element has been used for specifying the static forces. Following types of static 
forces are acting on the frame for all the three cases: 

Engine and gear box weight: 

Total weight of the engine and gear box is 102 Kgs. which is equally distributed on six nodes (six 
supporting points) of the shell elements on which engine and gear box is resting. 

Passenger load: 

Passenger load of 700 Kgs. is assumed on the platform for doing the analysis. On each node of 
the platform certain load has been assigned proportional to the area supported by it (calculated 
after generating the proper mesh). 
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pfq _ ”3^0 Freebodv diagram of sidemounting__model. 
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Gravity force: 


NISA provides facility to specify the body forces for the non-zero density structure. Gravity 
forces has been specified by using this facility for all the three cases. 


Dynamic forces: 

Engine unbalance forces: 

Engine unbalance forces as described in the Eqns. (3.4), (3.5) and (3.6) are applied equally on 
six nodes of the shell elements on which engine and gear hox is resting. Calculation of engine 
unbalance forces is done at following frequency: 

Idling: 

CASE I 

Idling frequency of the engine is 14.25 Hz. AH the unbalance forces for case I are calculated on 
this frequency. 

Ditch: 

CASE II 

For this case engine speed for calculating the engine unbalance forces has been determined as 
follows: 




A. 


51 



Distance traveled by the rear wheel when the front wheel 
passes over the ditch = v t 

= 11.11*0.09 
= 0.9999 m 
Wheel periphery = 7t D 

= 7t * 488 
= 1533.1 mm 

Number of revolution taken by rear wheel in 0.18 sec. = (2*999.9/1533.1) 

= 1.305 rev. 

Wheel speed = 7.247 revisec. 

Engine speed = 7.247*4.1*1.29 
= 38.3 Hz 

.-.Unbalance forces for case n are calculated for engine speed of 38.3 Hz. 

CASE ni 

In this case engine speed for calculating the engine unbalance forces has been determined as 
foUows: 

Length of the ditch curve = J ds 





1*0 09 

Jo 


(150co coscoi)^ + 
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By applying trapezoidal rule for time step of 0.045 sec. 
length of the ditch curve = 1053.29 mm 
Wheel periphery = it D 

= 71 * 488 
= 1533.1 mm 

Number of revolutions taken by rear wheel in 0.18 sec. = (2* 1053.29/1533.1) 

= 1.37 rev. 


.*. Wheel speed = 7.6288 rev ./sec. 

.'. Engine speed = 7.6288*4.1*1.29 
= 40.35 Hz 

Unbalance forces for case n and case IH are calculated for engine speed of 40.35 Hz. 


3.7 DISCREATIZATION OF THE CHASSIS FRAME AND DATA INPUT TO NISA 


For the finite element analysis the chassis is discretized into 254 elements and 280 nodes in the 
Side mounting model and into 236 elements and 241 nodes in the Central mounting model. The 
mass of the frame is applied through body force option of NISA. The various data input to NISA 
for the Side mounting frame is given in the tables 3.1 to 3.5 and for the Central mounting frame is 
given in the tables 3.6 to 3.10. Different input sections for modeling Side mounting frame are 
given in the Fig. 3.10 and for modeling Central mounting frame are given in the Fig. 3.11. 


53 



Prop. 

id. 


A 

(mm^) 


lyy 

(min'*) 


(mm'*) 


lyz 

(mm'*) 


TSF 

(mm) 


J 

(mm'*) 


Pt. no. 


SBFY 


SBFZ 


1 


477 


5.0e5 


7.29e4 


0.0 


1431 


841 


6.16e5 


6.16e5 


0.0 


3.1e-5 


1.23e6 


0.0 


0.0 


0.0 


0.0 


0.0 


4.7e-3 


0.0 


4.7e-3 


0.0 


0.0 


0.0 


0.0 


4.7e-3 


0.0 


4.7e-3 


0.0 


2827 


1.27e6 


1.27e6 


0.0 


6.27e5 


4 


0.0 


2.3e-4 


0.0 


2.3e-4 


2.3e-4 


0.0 


2.3e-4 


0.0 


1963 


6.13e5 


6.13e5 


0.0 


3.82e6 


0.0 


3.3e-4 


0.0 


3.3e-4 


3.3e-4 


0.0 


3.3e-4 


0.0 


490 


1.88e5 


1.88e5 


0.0 


l.Oe-4 


3.82e4 


0.0 


4.0e-3 


0.0 


4.0e-3 


4.0e-3 


0.0 


4.0e-3 


0.0 


200 


3.21e4 


2.78e4 


0.0 


200 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


21 


504 


4.5e5 


2.16e5 


0.0 


6.7e-5 


4.45 e5 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Table no. 3.1 Input data for element with NKTP = 12 (Side mounting) 


P.I. 

Ai (mm^) 

A 2 (mm^) 

6 1 

125 

125 


Table no. 3.2 Input data for element with NKTP = 14 (Side mounting) 
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P.I. 


7 

7.66e3 

8 

6.25e4 


Table no. 3.3 Input data for element with NKTP = 17 (Side mounting) 


P.I. 

Node no. 

t (mm) 

10 

1 

8 

10 

2 

8 

10 

3 

8 

10 

4 

8 


Table no. 3.4 Input data for element with NKTP = 20 (Side mounting) 


P.I. 

m (Kg.) 

11 

17 

12 

0.989 

13 

40.7 

14 

53.46 

15 

26.732 

16 

12.98 

17 

23.9835 

18 

69.877 

19 

91.787 

20 

45.89 


Table no. 3.5 Input data for element with NKTP = 26 (Side mounting) 
























Prop. 

id. 

A 

(mm^) 

lyy 

(inin‘‘) 


lyz 

(mm"*) 

TSF 

(mm) 

J 

(mm"*) 

Pt. no. 

SB¥Y 

SBFZ 


477 

5.0e5 





1 

0.0 

0.0 

6 

7.29e4 

0.0 

3 

1431 

2 

0.0 

0.0 








3 

0.0 

0.0 








4 

0.0 

0.0 


841 

6.16e5 





1 

0.0 

4.7e-3 

1 

6.16e5 

0.0 

3.1e-5 

1.23e6 

2 

4.7e-3 

0.0 








3 

0.0 

4.7e-3 








4 

4.7e-3 

0.0 








1 

0.0 

2.3e-4 

7 

2827 

1.27e6 

1.27e6 

0.0 


6.27e5 

2 

2.3e-4 

0.0 








3 

0.0 

2.3e-4 








4 

2.3e-4 

0.0 








1 

0.0 

3.3e-4 

2 

1963 

6.13e5 

6.13e5 

0.0 

- 

3.82e6 

2 

3.3e-4 

"o.o 








3 

0.0 

3.3e-4 


1 






4 

3.3e-4 

0.0 








1 

0.0 

4.0e-3 

3 

490 

1.88e5 

1.88e5 

0.0 

l.Oe-4 

3.82e4 

2 

4.0e-3 

0.0 








3 

0.0 

4.0e-3 








4 

4.0e-3 

0.0 








1 

0.0 

0.0 

10 

200 

3.21e4 

2.78e4 

0.0 

2 

200 

2 

0.0 

0.0 








3 

0.0 

0.0 








4 

0.0 

0.0 








1 

0.0 

0.0 

25 

504 

4.5e5 

2.16e5 

0.0 

6.7e-5 

4.45e5 

2 

0.0 

0.0 








3 

0.0 

0.0 








4 

0.0 

0.0 








1 

0.0 

3.6e-3 

24 

575 

2.68e5 

2.68e5 

0.0 

5.7e-5 

5.34e5 

2 

3.6e-3 

0.0 








3 

0.0 

3.6e-3 








4 

3.6e-3 

0.0 


Table no. 3.6 Input data for element with NKTP = 12 (Central mounting) 


P.I. 

Ai (mm^) 

Az (mm^) 

4 

125 

125 


Table no. 3.7 Input data for element with NKTP = 14 (Central mounting) 
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P.I. 

K(k 2 /s") 

5 

7.66e3 

8 

6.25e4 


Table no. 3.8 Input data for element with NKTP = 17 (Central mounting) 


P.I. 

Node no. 

t (mm) 

9 

1 

8 

9 

2 

8 

9 

3 

8 

9 

4 

8 


Table no. 3.9 Input data for element with NKTP = 20 (Central mounting) 


P.I. 

ra (Kg.) 

11 

17 

12 

4.1 

13 

10.05 

14 

11.94 

15 

13.125 

16 

14.31 

17 

7.2 

18 

40.55 

19 

38.16 

20 

30.29 

21 

28-5 

22 

70.805 

23 

66.66 


Table no. 3.10 Input data for element with NKTP = 26 (Central mounting) 
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Fig. 3.10 Different input sections for sidemounting frame 







CHAPTER - 4 


THEORETICAL RESULTS 


4.1 MAXIMUM SHEARING STRESS THEORY 


Observations made in the course of extrusion tests on the flow of soft metals through orifices lend 
support to the assumption that the plastic state in such metal is created when the maximum 
shearing stress just reaches the value of the resistance of the metal against shear. Let cti , 02 and as 
be the three principal stresses acting on three mutually perpendicular planes at a point in an 
isotropic material. Assuming ai>a2>a3 , yielding , according to this theory occurs when the 
maximum shearing stress 


X 


max 




2 


reaches a cntical value. The maximum shearing stress theory is accepted to be fairly justified for 
ductile materials i.e. this theory says that if ai>a2>a3 are the three principal stress at a point, 
failure occurs when 


0^1 




2 


> 



2 


where ay /2 is the shear stress at yield point m a uniaxial test 
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4.2 RESULTS 


NISA output gives us the maximum sheanng stresses for maximum of four points on a given 
cross section. Since it is a 3-D combined bending and torsion problem so the maximum stress for 
a NISA beam element will occur on the surface of the element By doing preliminary analysis 
four points on each cross section have been identified (Fig. 3.10 and Fig. 3.1 1) to represent that 
cross section for analysis for critical stresses. For Case I stress analysis has been done for one 
engine cycle and for Case II and Case IH stress analysis has been done for the time for which 
vehicle passes over the ditch. In all the three cases total time span has been divided into eight 
steps for software to give the output. Alter sorting out the maximum sheanng stresses on each 
cross section for selected nodes table has been formed for the sorted data. A computer program 
has been written to get the time Vs factor of safety data using maximum shearing stress theory. 
This data has been plotted using a software package. 

Since it is a general 3-D structural dynamic force analysis it is expected that the maximum stress 
point at cross section will vary with time. This is verified by the data in the tables 4.1 to 4.6. At 
each cross section there is no fixed maximum stress point. Maximum stress point varies with time 
for a cross section. Also the time at which the maximum stress on a cross section for the whole 
duration of analysis occurs is different for different cross sections. 

For Case I when the engine is idling and vehicle is stationary, maximum shearing stress for Side 
mounting frame occurs at node 273 for section 1, at node 59 for section 2, at node 14 for section 


63 



29 

^ O ; 

^ Cl 

JO 

S S c 

^ < 
N G ; 

^ d < 

o 

■M 

S d 

o 

g G ; 

■ 

34 

xr 

ro 

cs 

voj Cl 
»«— -< 


ro 

^ G 

CO 

CO 

G 

04 

r4 

o G 

cr\ 

CO „ 

c^ d ■ 

r— < 



d 

04 

VO 

CO 

o ^ 

o ^ 

CO Cl 

CO 3^ 

3 G 

CO ^ 

^ G 

r-H 

04 

d 

04 d 

O' 

CO ^ 

5 G 

xO ^ 

<0 04 

CO ' 



(0 "TT** 

04 Cl 

04 

CO _ 

xo 31 

CO 

CO ^ 

- G 

CO ^ 

G 

04 

c^ ^ 

<2^ dr* 

CO 31 

CO ^ 

VO G 

O^ 

CO ^ 

^ G 

CO 

270 

CO ^ 

oi 

CO _ 

O' Cl- 
vn 

CO ^ 

dr 

c\ cl 

to 

CO _ 

o ^ 

CO 

04 _ 

^ G 

O- 

CO ^ 

CN G 

T i 

CO ^ 

VO G 
■*— -< 

CO ^ 

O CO 

'::t — 

59 

CO 

0 ^ 

01 

CO „ 

O 

O- 

vo> 

CO 

<0 t--< 

VO 

CO ^ 

dr 

C4 3l 
CO 

04 _ 

(D 

CO 31 
r- 

CO ^ 

c^ d 

CO _ 

2 G 

CO ^ 

d 

CO 

CO 

Xl- 

o ^ 

XO ' 

o4 - 

^ S 
<0 

CN _ 

<o ^ 
to 

VO 

C4 _ 

^ dr 

o d 

04 _ 

o :5^ 
VO 31 

XT 

C4 _ 

G 

O' 

04 

G 

CO 

04 ^ 

CO 

45 

04 ^ 

xr 

04 

VO G- 
xf 

04 ^ 

^ d 

04 _ 

^ 04 

c\ d 

C4 _ ■ 

d 

to 

C^4 

^ d 

VO 

04 _ 

d 

CO 

04 ^ 

c^ G 

VO 

273 

04 

04 ^ 

O — ^ 

d 

CO ^ 

^ d 

CO ^ 

^ d 

t~-< 

CO 

s s 

04 ^ 

^ CO 

CO d 

C\ 

04 ^ 

G 

r 

04 ^ 

C§^ G 

xO 

■oil 

J.755C-2 

04 

1 

o 

u-> 

CO 

04 

6 

to 

VO 

<N 

to 

C4 

t 

o 

C4 

o 

C-4 

1 

O 

to 

r- 

cd 

T— « 

1 

o 

CO 

to 

o 

1 

o 

CO 

04 

04 

1 

o 

- ^ 
o 

XT 


64 


Tabic no. 4,1 Varialion of max, shearing stress (Kg/ s’/mm) with lime - Case I (Side mountmg) 
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3.1 and at node 34 for section 3.2. For the same case maximum shearing stress for the Central 
mounting frame occurs at node 34 for section 1, at nodes 60 and node 56 for sections 2.1 and 2.2 
respectively, at nodes 83 and 81 for sections 3.1 and 3.2 respectively, and at node 78 and node 
74 for sections 3.3 and 3.4 respectively. Except for section 2 maximum shearing stress occurs at 
different points in space in both the frames for sunilar sections. Value of maximum shearing stress 
is quiet high for sections 3.3 and 3.4 in the Central mounting frame which makes it a weak point 
for failure. The sunilar section is not present in the Side mounting frame. Fig. 4.2 also shows that 
factor of safety remains very low for sections 3.1 to 3.4 of Central mounting frame for most of the 
loading duration. 

For Case n when the vehicle is running at 40 kmThr. in the fourth gear and its front wheel goes 
into a ditch of half sine wave of 2 m wavelength and 150 mm amplitude, maximum shearing stress 
for Side mounting frame occurs at node 49 for section 1, at node 60 for section 2, at node 14 for 
section 3.1 and at node 34 for section 3.2. For the same condition maximum shearing stress for 
the central mounting frame occurs at node 34 for section 1, at nodes 60 and 56 for sections 2.1 
and 2.2 respectively, at nodes 83 and 81 for sections 3.1 and 3.2 respectively, and at nodes 78 and 
74 for sections 3.3 and 3.4 respectively. In this case also the location of the maximum shearing 
stress points are same in space as is in the previous case, for section 2 and section 3 of both the 
models. Fig. 4.4 once again shows that section 3.3 and section 3.4 have lower factor of safety for 
most of the loading duration and are poor sections. 
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For Case DI when the vehicle is running at 40 km./hr. in the fourth gear and its both rear wheels 
go into a ditch of half sine wave of 2 ni wavelength and 150 mm amplitude, maximum shearing 
stresses for the Side mounting frame occurs at node 274 for section 1, at node 60 for section 2, 
at node 14 for section 3.1 and at node 34 for section 3.2. For the same conditions for the Central 
mounting frame maximum shearing stresses occurs at node 28 for section 1 , at nodes 60 and node 
56 for sections 2.1 and 2.2 respectively, at nodes 78 and 74 for sections 3.3 and 3.4 respectively. 
Once again in this case also for section 2 and section 3 maximum sheanng stress points are 
located at the same place in space as is in the previous two cases, for both the two models. Again 
sections 3. 1 to 3.4 of the Central mounUng frame have the larger maximum shearing stress and 
sections 3.3 and 3.4 have the largest maximum shearing stress values. 

From the above results it can be concluded that in both the models section 1 has the lowest stress 
developed. Section 1 and section 2 in both the models develop similar stress values for aU the 
three cases. But section 3 in the Central mounting frame invariably has higher stress values in 
comparison to section 3 for the Side mounting frame. Sections 3.3 and 3.4 of the Central 
mounting frame are the weakest sections . Analysis shows that inspite of the Central mounting 
being a more balanced frame and engine being mounted in the middle of the frame , is a weaker 
frame due to sections 3.3 and 3.4. Similar sections are not present in the Side mounting frame. 
This can be attributed to the fact that the diameter of the load carrying tube has been reduced and 
the channel section has been replaced by tubural section in the Central mounting frame. This 
change reduces the area moment of inertia of load carrying member and hence increases the 
stresses developed in it. 
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CHAPTER - 5 


EXPERIMENTAL VERIFICATION 


5.1 THEORY OF STRAIN GAGE 


The electrical resistance strain gage is basically a piece of very thin foil of fine wire Fig. 5.1 which 
exhibits a change in resistance proportional to the mechanical strain imposed on it. In order to 
handle such a delicate filament, it is either mounted or bonded to some type of carrier material and 
is known as the bonded strain gage. The strain gage is used universally in the experimental 
determination of stresses. The relation ship between resistance change and strain in the foil or 
wire used in strain gage construction can be expressed as: 


AL _ AR 
L ~ k*R 


or 
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Fig. 5.1 Typical construction of a foil strain gage 
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whdre k is defined as the gage factor of the foil or wire, AR is the resistance change due to strain , 
R is the initial resistance, AL is the change in length, L is the original length of the wire or foil, 
and AL/L is the unit strain to which the wire or foil is subjected. 

5.2 DETAILS OF BRIDGE CIRCUIT 

The bridge circuit converts the change in resistance of the strain gauges bonded to the chassis 
frame into a potential difference. This potential difference is in turn, recorded on the oscilloscope. 
Fig. 5.2 shows the configuration of the bridge circuit used in this smdy. Ri ,R 2 , R 3 and R 4 are 
four strain gauges each having a resistance value of 120 ±0.3 £2 . Ri represent the strain gauges 
bonded to the chassis frame. R 2 ,R 3 and R 4 are the dummy gauges. These dummy strain gauges 
are bonded on a steel plate which acts as a heat sink to the heat generated in the cods of the strain 
gauge when current flows through it. The bridge circuit is balanced to give zero voltage output 
for no change in resistance of the active strain gauges. This is done by connecting a 1.0 £2 
resistance in series and two variable resistors (1.0 M£2 and 1.0 k£2 ) parallel to one of the 
dummy strain gauges (R 2 in the shown circuit). To calibrate the bridge circuit a calibration 
resistance R^ is connected parallel to the active gauge through a switch ( K ). The value of the 
calibration resistance is 47.0 k£2. The equation involving the calibration of the bridge circuit are as 
follows; 

For a balanced bridge circuit the output voltage ( Ae ) is given by: 
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where E is the input voltage. It can be seen from the equation that similar ( both positive or both 
negative ) changes in resistance of opposite arms of the bridge circuit are added up and dissimilar 
( one positive and other negative ) changes are canceled out. 

The relation between the strain in the strain gauge and corresponding change in its resistance is 
governed by the Eqn. (5.1). Which is: 


G.F. 

L 

The change in resistance of the arm CD after connecting Rc is given by 


AL/ 


AR/ 




G.F. 


AR = R- 


R + R^^ 


leading to 



R 1 : Active gage 

R R ^ , R ^ : Dummy gages 



AR _ R 
R R + R 


(5.3) 


corresponding to this change in resistance a voltage difference ( calibration voltage, Vc ) will 
occur between terminals A and C. This resistance change will also correspond to the strain in the 
strain gauge. From Eqns. (5.1 ) and (5.3) this strain will be: 


_ AL_ R 
" L ~ iR + R^)*G.F. 


(5.4) 


The strain value £c will thus correspond to the calibration voltage, Vc.. By the linearity of the 
relation ( Eqns. (4.1) and (4.2) ) between voltage drop across AC strain in the strain gauges 
bonded to the frame, the strain in the frame corresponding to a voltage V recorded on the 
oscilloscope can be given as: 


e = 



(5.5) 


5.3 OSCILLOSCOPE 


1604 Digital storage oscilloscope has been used for experimental work. It can be operated in 
storage or non-storage modes as required. Two complete traces can be stored for future use and 
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four complete front panel control setup. With the optional battery backup facility, these can be 
retained even when it is switched off. Its input sensitivity is 500jJ.V to lOV/div. This osciUoscope 
is provided with differential input facihty. It has a band width of 3 MHz which means it can record 
signals up to 3 MHz frequency. 

5.4) DETERMINATION OF PRINCIPAL STRAIN AND ITS DIRECTION: 


For doing the verification of the models developed for NISA for doing analytical analysis, 
experiments have been performed on the chassis. Strain gauge has been put in the maximum 
absolute principal stress direction to measure the strain . This experimental value of strain has 
been compared with the theoretical value of strain given by theoretical analysis. Theoretically 
strain has been calculated from the following formulae: 

( 5 . 6 ) 

E ^ 

Where, 

El : Principal strain at a point 

CTi , ( 72 , as : Three principal stresses acting on three mutually perpendicular planes at a point 
v: Poisson’s ratio 
E: Young’s modulus 

For determining the principal strain direction following theory has been developed: 
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From 3-D beam bending theory stress tensor is given as: 


0 0 

5 ^' 0 0 


and from torsional theory stress tensor is given as: 


0 V 
V 0 0 

0 0 


(5.7) 


(5-8) 


3-D beam element in NISA incorporates both the bending and the torsion. Since element of the 
stress tensors can be added up to combine the effect of bending and torsion hence for 3-D beam 
element after adding Eqn.(5.7) and Eqn.(5.8) we get the stress tensor as: 


[a] = 



0 0 
0 0 


(5.9) 


NISA gives all the elements of Eqn.(5.9), principal stresses and maximum shear stress in the 
output. Since the chassis is made of isotropic material (no directional property) the direction of 
the principal strains and principal stress will coincide. To calculate principal strain direction let us 
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assume that there is a plane n with direction cosines n* , Uy and on which the stress is wholly 
normal. Then we have: 


T = a n 


The components of this along the x, y and z axes are: 


Tx = <7 nx , Ty = a ny and Tz = a nz 


Stress tensor Eqn.(5.9) can also be written as: 


Sxx nx + Sxy Uy + Sxz nz — Tx 


Sxy Rx — Ty 


Sxz Rx — Tz 


from Eqn. (5.11) and Eqn. (5.12) we get: 



(5.10) 

(5-11) 

(5.12.1) 

(5.12.2) 

(5.12.3) 

(5.13.1) 


86 



(5.13.2) 



We also know that 


2 2 “2 

Hx + ny^ + Hz 


(5-14) 


on solving Eqn. (5.12) and Eqn. (5.13) 


^ 2 ^ 1 

l + (^)^ + (^)2 
a a 




1 


•S'xv 7 *5 , 

l4.(_5L)2 + (fs.)2 

a G 


(5.15.1) 




“'xy 


' G G 


(5.15.2) 


ajl + (^f + (^y 

V cr G 


(5.15.3) 


5.5 EXPERIMENTAL SETUP 


Fig. 5.3 shows the schematic diagram and Fig. 5.5 shows the instrumentation for the complete 
experimental setup. Strain gage is pasted, Fig. 5.4, on the chassis at the critical points in the 
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Fig. 5.3 Schematic diagram of experimental setup 
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Fig. 5.4 Strain gage pasted on chassis member 


Fig. 5.5 Oscilloscope with bridge circuit box 
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principal strain direction. It is connected to the bridge circuit through long coaxial wire. Bridge 
circuit gives the output to the oscilloscope. 


5.6 RESULTS 


Experiments has been performed for Case I i.e. for the condition when engine is idling. Table 5.1 
and table 5.2 shows the expected theoretical maximum principal strain and the experimental 
principal strain values for selected nodes . for which experiments has been performed. 

Maximum error in the expected value and experimental value is 94.7 %. Error varies from 8 % 
to 94.7 %. This error can be attributed to the simplifications made in the computer modeling due , 
to various reasons e.g. inability of the package to model the engine, power transmission system, 
tyres, non availability of the data for unbalanced engine forces, shock absorbers, complications 
involved in modeling for welded plates, distribution of load etc.. Error can also be attributed to 
the fact that experimental values has been taken when the steady state condition has been reached 
where as for doing computer modeling because of non availability of the steady state boundary 
condition data simphfied boundary conditions have been applied. 


Since various simplifications are done to do the computer modeling it is expected that the 


theoretical strain values will be more than the experimental strain values, which is \erified by the 


principal strain values in the table 5. 1 and table 5.2. With the limitation imposed for doing the 
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Node no. 

Theoretical 

strain value 

Experimental 

strain value 

% Error 

270 

9.79e-4 

1.989e-4 

79.7 

14 

3.44e-4 

0.37e-4 

. 89 


Table no. 5.1 Comparison of principal strains (Side mounting) 


Node no. 

Theoretical 

Experimental 

% Error 


strain value 

strain value 


62 

2.46e-4 

0.54e-4 

78 

65 

8.68e-4 

0.46e-4 

94.7 

78 

11.5e-4 

10.625e-4 

8 


Table no. 5.2 Comparison of principal strains (Central mounting) 
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computer modeling it can be said that the errors are with in the acceptable norms and the 
computer model is fairly good for doing the required analysis for getting the qualitative solution. 
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CHAPTER - 6 


CONCLUSION AND SCOPE FOR FUTURE WORK 


In the present work computer modeling of the two chassis (with old and new design) acted upon 
by d 3 mamic and static loads has been done. Finite element analysis of the two chassis models is 
done by using these computer models, to determine and compare stresses at critical sections under 
the loading conditions. These finite element analysis model and results have been verified by 
expenmental stress analysis. 

In this work various simplifications have been done due to various reasons. By using abetter 
software package these simplifications can be reduced and stress predications can be done with 
greater accuracy. In the present work finer mesh can be generated and more points can be 
selected for finer stress variation analysis. 


For doing experimental analysis only a few points have been selected. The number of points for 
doing experimental analysis can be increased. Experimental stress values can also be taken at 
different engine speeds. By usmg strain rosettes strain fields at various critical points can be found 
out and detailed experimental stress analysis can be done. For Case 11 and Case HI theoretical 
results can be verified by using a pulsator to simulate the ditch boundary condition. Also vehicle 
can be run in the actual road conditions and by using on board instruments experiments can be 
performed for real conditions. 
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APPENDIX 


dimension a(8),b(8) 

a( I ) = 0 . 01 1 25 
b( 1) = 0.01755 
sum! = 0.01125 
sum2 = 0 . 01755 
do i = 1,7 
a( i+1 ) = aC i ) + sual 
b( i+1 ) “ bC i ) + suaS 
enddo 

opentunit « 10, file = 
open(unit = £0,file = 

X = £20000 

do i = 1,8 
read ( 1 0 , ♦ ) y 
s = x/y 

wriTe(20,»)a(i),s 

wriTe{£0,*)b(l),s 

enddo 

stop 

end 


'cd81 .in' ) 
'NodeSl ' ) 


97 



